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Abstract 

We report a measurement of the ratios of Wilson coefficients A9/A7 and Aiq/Aj in B — > K*l + £~ . 
The result is obtained from a data sample containing 386 million BB pairs that was collected at 
the T(4S I ) resonance with the Belle detector at the KEKB asymmetric energy e + e" collider. 

PACS numbers: 13.20.He, 11.30.Hv 
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Flavor-changing neutral current (FCNC) b — > s processes are forbidden at tree level in 
the Standard Model (SM); rather, they proceed at a low rate via loop or box diagrams. 
If additional diagrams with non-SM particles contribute to such a decay, their amplitudes 
will interfere with the SM amplitudes and thereby modify the decay rate as well as other 
properties. This feature makes FCNC processes an ideal place to search for new physics. 

To evaluate the new physics contributions in b — > s processes, Wilson coefficients are 
used 0. These coefficients parameterize the strength of the short distance interaction. If 
new physics contributes to the b — > s processes, the relevant coefficients will deviate from the 
SM values. For electroweak penguin decays, the effective Wilson coefficients Cy 5 , Cg S and 
C^q appear in the partial decay width. A next-to-next-to-leading order (NNLO) calculation 
for these effective coefficients has many correction terms [2], so leading coefficients A 7 , Ag 
and A w are usually used for the evaluation. 

Measurements of the radiative penguin decay B — > X s ^ 0, 0, H|, which are consistent 
with the SM prediction, strongly constrain the magnitude — but not the sign — of the Wilson 
coefficient A 7 . This is an important constraint; however, non-SM contributions can change 
the sign of A 7 without changing the B — ► X s ^ branching fraction |(|. 

The b — > s£ + £~ process is promising from this point of view, since not only the photonic 
penguin diagram but also the Z-penguin and box diagrams contribute to this decay mode. 
As a result, we can determine the relative signs of the Wilson coefficients A 7 , Ag and Axo 
as well as their absolute values. The first observations of B — > K£ + £~, B — > K*£ + £~ 
and inclusive B — > X s £ + £~ decays were reported by the Belle Collaboration 0, H, 0| • The 
measured branching fractions of these decay modes E , pi Ijj. fiol 11 1 were used to exclude a 
large area of the allowed region in the Ag-A w plane [l2 . fl3Ul4j] . However, the determination 
of the relative sign of A 7 to Ag or A w (as well as Ag and A w ) requires precise measurements 
of differential branching fraction and the forward-backward asymmetry as functions of q 2 
in these decay modes [15| . As shown in Fig. [H the distributions are quite sensitive to the 
choice of Wilson coefficients. 

The first measurement of the forward-backward asymmetry in B — ► K^£ + £~ was re- 
ported by the Belle Collaboration PH. The result for B -> K*£+£~ is consistent with both 
the SM case and the case of a flipped sign for A 7 Ai , however a quantitative evaluation in 
terms of Wilson coefficients has not been done yet. 

In this paper, we present preliminary results of a measurement of Wilson coefficients in 
B — > K*£ + £~ using data produced in e + e~ annihilation at the KEKB asymmetric collider 
(l7| and collected with the Belle detector. The data sample corresponds to 357 fb _1 taken 
at the T(49) resonance and contains approximately 386 million BB pairs. As a control 
sample, we also study the B + — > K + £ + £~ mode which has very small forward-backward 
asymmetry regardless of the existence of new physics Here and throughout the paper, 
the inclusion of charge conjugate modes is implied. 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cerenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter comprised of CsI(Tl) crystals (ECL) located 
inside a super-conducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside of the coil is instrumented to detect K\ mesons and to identify muons 
(KLM). The detector is described in detail elsewhere jl^ . Two inner detector configurations 
were used. A 2.0 cm beampipe and a 3-layer silicon vertex detector was used for the first 
140 fb _1 data sample (DS-I) containing 152 million BB pairs, while a 1.5 cm beampipe, a 4- 
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FIG. 1: Expected forward-backward asymmetry and q 2 distributions for several input parameters. 
The discontinuity at 5.5 GeV/c 2 is due to the imperfect matching between full and partial NNLO 
calculations. 

layer silicon detector and a small-cell inner drift chamber were used to record the remaining 
217 fb~ x data sample (DS-II) containing 234 million BB pairs [2oj | . 

In this analysis, primary charged tracks — except for the daughters from Kg — > 7r + 7r~ 
decays — are required to have impact parameters relative to the interaction point of less 
than 5.0 cm along the z axis (aligned opposite the positron beam) and less than 0.5 cm 
in the rip plane that is transverse to this axis. This requirement reduces the combinatorial 
background from photon conversion, beam-gas and beam-wall events. 

Charged tracks are identified as either kaons or pions by a likelihood ratio based on the 
CDC specific ionization, time-of-flight information and the light yield in the ACC. This 
classification is superseded for a track that is identified as an electron or for a pion-like track 
that is identified as a muon. Electrons are identified from the ratio of shower energy of the 
matching ECL cluster to the momentum measured by the CDC, the transverse shower shape 
of this cluster, the specific ionization in the CDC and the ACC response. We require that 
the electron momentum be greater than 0.4 GeV/c to reach the ECL. Muons are identified 
by their penetration depth and transverse scattering in the KLM. The muon momentum 
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is required to exceed 0.7 GeV/c. The muon identification criteria are more stringent for 
momenta below 1.0 GeV/c to suppress misidentified hadrons. We also apply a kaon veto for 
electron and muon candidates. 

Photons are selected from isolated neutral clusters in the ECL with energy greater than 
50 MeV and a shape that is consistent with an electromagnetic shower. Neutral pion can- 
didates are reconstructed from pairs of photons, and are required to have an invariant mass 
within 10 MeV/c 2 of the nominal 7r° mass and a laboratory momentum greater than 0.1 
GeV/c. The pion momentum is recalculated by constraining the invariant mass to the nom- 
inal 7T° mass. Kg candidates are reconstructed from oppositely charged pions that have 
invariant masses within 15 MeV/c 2 of the nominal Kg mass. We impose additional criteria 
based on the radial impact parameters of the pions (5r), the distance of closest approach of 
the pions along the beam direction (Sz), the distance of the vertex from the interaction point 
(I), and the azimuthal angle difference between the vertex direction and the Kg momentum 
direction (5(fi). These variables are combined as follows: for p K o < 0.5 GeV/c, we require 
5z < 8 mm, 5r > 0.5 mm, and 5(fi < 0.3 rad; for 0.5 GeV/c < p K o < 1.5 GeV/c, we apply 
5z < 18 mm, 5r > 0.3 mm, 5(f) < 0.1 rad, and I > 0.8 mm; and for p K o > 1.5 GeV/c, we 
impose 5z < 24 mm, 8r > 0.2 mm, 5<p < 0.03 rad, and I > 2.2 mm. 

K* candidates are formed by combining a kaon and a pion: K + tt~, K s tt + or K + tt°. 
The K* invariant mass is required to lie within 75 MeV/c 2 of the nominal K* mass. For 
modes involving neutral pions, combinatorial backgrounds are reduced by the additional 
requirement cos^hei < 0.8, where 6*hei is defined as the angle between the opposite of the B 
momentum and the kaon momentum direction in the K* rest frame. 

B candidates are reconstructed from a K^*> candidate and an oppositely-charged lepton 
pair. If a photon with energy less than 500 MeV is found in the 50 mrad cone along the 
electron or positron momentum direction, we add the four momentum of the bremsstrahlung 
photon to the di-electron system. We use two variables defined in the CM frame to select 



B candidates: the beam-energy constrained mass M bc = J E^ eSLm 2 — p* B 2 and the energy 
difference AE = E* B — E^ eam , where p* B and E* B are the measured momentum and energy, 
respectively, of the B candidate, and E^ e&m is the beam energy |21j. When multiple candi- 
dates are found in an event, we select the candidate with the smallest value of |Ai?|. 

Backgrounds from B — > J/ip(ip')K^ are rejected using the dilepton invariant mass. The 
veto windows are defined as 

-0.25 GeV/c 2 < M ee - M J/V , < 0.07 GeV/c 2 for K* modes 

-0.20 GeV/c 2 < M ee - M m < 0.07 GeV/c 2 for K modes 
-0.20 GeV/c 2 < M ee - My < 0.07 GeV/c 2 for K* and K modes 
-0.15 GeV/c 2 < M w - M J/V , < 0.08 GeV/c 2 for K* modes 
-0.10 GeV/c 2 < M w - Mj/^ < 0.08 GeV/c 2 for K modes 
-0.10 GeV/c 2 < M w - My < 0.08 GeV/c 2 for K* and K modes 
If a bremsstrahlung photon is found, we reapply the above vetoes with the invariant mass cal- 
culated including this photon, to reject the background J/i/j(i/j') —>■ e + e~ 7(7). For K*e + e~ 
modes, B — ► J /ip K can be a background if a bremsstrahlung photon is missed and a pion 
from the other B meson in the event is included. We suppress this background by the 
following prescription: the included pion is discarded and an unobserved bremsstrahlung 
photon is added with a direction parallel to the electron or positron and an energy that 
gives AE = for the B candidate. If the dilepton mass and the beam-energy constrained 
mass are consistent with a B — > J/ipK event, the candidate is vetoed. 

We suppress background from photon conversions and 7r° Dalitz decays by requiring the 
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di-electron mass to satisfy M e + e - > 0.14 GeV/c 2 . This cut eliminates possible peaking 
background from B — > K*^ and B — > K^n°. 

Background from continuum qq events is suppressed using event topology. Continuum 
events have a jet-like shape while BB events have a spherical shape in the center-of-mass 
frame. A Fisher discriminant T jl^ is calculated from the energy flow in 9 cones along the 
B candidate sphericity axis and the normalized second Fox- Wolfram moment R2 jlH • We 
combine this with the cosine of the polar angle cos 9* B of the B meson flight direction and the 
cosine of the polar angle cos * h of the B meson sphericity axis to define likelihoods £ S ig and 
Axmt for signal and continuum background, respectively, and then cut on the likelihood ratio 
= £ s i g /(£ s i g + £ cont ). F° r the muon mode, | cos6>* ph | is not used since its distribution 
is nearly the same for signal and continuum within the detector acceptance. 

The dominant background from BB events is due to semileptonic B decays. The missing 
energy of the event, E miss = 2E£ eam — E* is where E* is is a total visible energy in the event, 
is used to suppress this background since the undetected neutrinos carry away a substantial 
amount of energy. The B meson flight angle, cos8 B , is also used to suppress combinatorial 
background in BB events. We combine -Emiss and cos 6^ into signal and i?i?-background 
likelihoods and cut on the likelihood ratio TZ BB , defined similarly to 1Z cont . 

The signal box is defined as |Mb c — m B \ < 0.008 GeV/c 2 for both lepton modes and 
-0.055 GeV< AE < 0.035 GeV (\AE\ < 0.035 GeV) for the electron (muon) mode. We 
make distinct selections on TZ CO nt and 7Z BB - for each decay mode to optimize the sensitivity 
for the region with q 2 less than 6 GeV/c 2 . 

To determine the signal yield, we perform an unbinned maximum-likelihood fit to the 
M bc distribution of events with M bc > 5.2 GeV/c 2 and -0.055 GeV< AE < 0.035 GeV 
(|A£J| < 0.035 GeV) for the electron (muon) mode. We consider the signal, three cross-feed 
and four background components. The cross-feed components are correctly and incorrectly 
flavor tagged K*£ + £~ events ( "cfcf " and "ifcf " , respectively) and b —> s£ + £~ processes other 
than K*£ + £~ (X s £ + £~ events). The four background components are dilepton background, 
-fT*£ ± ^ =F background, K*h + h~ background and ip background, where h refers to a pion or a 
kaon. The dominant components in the signal box are signal and dilepton background. For 
the muon mode, K*^^ background is sub-dominant. The ratio of K*^b^ to dilepton 
background is about 0.24 and 0.12 for muon and electron modes, respectively. 

The expected number of signal events is calculated as a function of M^ c using a Gaussian 
signal distribution plus background functions. The mean and the width of the signal Gaus- 
sian are determined using observed J /ipK^ events. A MC study shows that the width has 
no dependence on the dilepton invariant mass. The shape and normalization of cross-feed 
events are parameterized by a Gaussian and an ARGUS function [24J based on MC samples. 

The background from semileptonic decays is parameterized by an ARGUS function. The 
shape of the backgrounds are determined from large BB and continuum MC samples in 
which each event contains at least one oppositely charged lepton pair. The shape parame- 
ter obtained from MC is consistent with that taken from a data sample of B — > K^e ± fx T 
candidates. Since the M^ c shape for K*£ background is similar to that for the dilepton 
background, we use the same shape function. The residual background from J/ip and if}' 
mesons that cannot be removed by the if>( ) veto windows is estimated from a large MC 
sample of J/if> and if)' inclusive events and parameterized by an ARGUS function and a 
Gaussian. The background contribution due to misidentification of hadrons as leptons is 
parameterized by another ARGUS function and a Gaussian. The ARGUS function repre- 
sents the combinatorial background while the Gaussian represents the component that forms 
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a peak in the M bc distribution. The shape and normalization of this background are fixed 
using the B — > K^h + h~ data sample (where ft, refers to a pion or kaon). All K^h + h- 
combinations are weighted by the momentum- and polar angle-dependent probability of 
misidentifying K^h + h~ as K^£ + £-. This study yields 3.8 ± 0.4 and 2.5 ± 0.3 peaking 
events in the signal box for Kh + h~ and K*h + h~ , respectively. Other backgrounds with 
misidentified leptons are negligible. The probability density function(PDF) for the M bc fit 
is expressed as follows: 

P(M hc ) = f sig G sig (M hc ) 

+ /sig/cfcf (fg { G cM (M hc ) + (1 - /£ f M cfcf (M bc )) 

+ /sig/ifcf (/£ f G ifcf (M bc ) + (1 - ^ cf )Af C f(M bc )) 

+ (fleeG XaU (M hc ) + (1 - /£«)A Xs «(M bc )) 

+ /dilep^dilep(^bc) 

+ fK'kh (fg* hh G K * hh (M bc ) + (1 - /g.jAjc.^CMbe)) 

+ U (/^(M bc ) + (1 - /«)^(M bc )) , (1) 

where / is the fraction of each of signal and backgrounds, G is a Gaussian and A is an 
ARGUS function. The / S j g and /dii ep are floated in the fit. 

Figure El shows the fit result. We obtain 96.0 ± 12.0 and 113.6 ± 13.0 signal events for 
K + £ + £~ and K*£ + £~, respectively. Table summarizes the signal yield for each subdecay 
mode. 



Mode 


Signal yield 


K+e + e- 


46.6 ±8.7 


(K + 7r-)*°e + e- 


22.4 ±7.0 


(K° s TT + )*+e + e- 


10.9 ±4.4 


(K+ir°)*+e+e- 


5.8 ±4.0 




49.4 ±8.2 




56.6 ±8.8 


(^7T+)*+/X+ M - 


7.7 ±4.0 




10.2 ±4.1 



TABLE I: Signal yields obtained from the M bc fit. The errors are statistical only. 

To extract Wilson coefficients, we perform an unbinned maximum likelihood fit to the 
normalized double differential decay width (l/r)<ir/(ig 2 c?cos6 l for B — > K*£ + £~ [25J, where 
cos 9 is the cosine of the angle between negative (positive) charged lepton and B° or B + (B° 
or B~) meson momenta in the dilepton rest frame. The forward-backward asymmetry is 
defined as 

2 r(g 2 ,cosg>0)-r(g 2 ,cosg<0) 
B r(g 2 ,cos^>0)±r(g 2 ,cos^<0)' [) 

We fix all parameters in the to be SM values except the leading coefficients A 7 , A 9 and 
A w . In the SM at renormalization scale /i = 2.5 GeV, the values of the leading coefficients are 
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FIG. 2: M bc distributions for (a) B -> if+e+e", (b) B -» K*e + e -, (c) B -» 2f+ /i+ p~ , (d) 
-B — > K*fi + {i~, (e) -B — > and (f) -B — > K*l + l~ samples. The solid and dashed curves are 

the fit results of the total and background contributions. 
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A 7 = -0.330, A 9 = 4.069 and A 1Q = -4.213 HQ. In the double differential decay width, 
there are a few additional parameters to be fixed. These include the form factor parameters 
and the bottom quark mass m&. We choose the form factor model of Ali et al. [l2l . 25[ and 



the bottom quark mass of 4.8 GeV/c 2 for the nominal fit. The renormalization scale \x is 
set at 2.5 GeV as suggested by Ref. [l2T |. 

The PDF used for the fit consists of terms describing signal, cross-feeds and backgrounds: 

P{q\ cos6, M bc ) = / sig (M bc )e? g V, ™s6) ^^ (q 2 , cos 9; A 7 , A 9 , A W )/N™(A 7 , A 9 , A w ) 

+ /cfcf(M bc )e cfcf (g 2 , cos ^ dq 2^ cos Q ^ q2 ' cos6l; A w )/N cM (A 7 , A 9 , A w ) 

d 2 T 

+ /ifcf(M bc )e ifcf (g 2 , cos6) — -(q 2 , - cos#; A 7 , A Q , A 10 )/N m {A 7 , A 9 , A 10 ) 

dq l d cos 

+ fxAM hc )V Xs ee(<l 2 ^os6) 

+ /diiep(M bc ){(l - fK*eh)Vdn cp {q 2 ,cos6) + f K *ehV K *eh{q 2 , cos 0)} 

+ fK*hh(M hc )V K * hh (q 2 , cos 6) 

+ U{M hc )V4q 2 ,cos6), (3) 

where / and V are the event by event fraction and probability density function for each 
components, respectively. N is the normalization factor and e is the efficinecy function. Each 
of the fractions (/ S i g etc.) are determined mode by mode from the M bc fits except fK*eh, 
the fraction of K^l^h^ background within the dilepton background component, which is 
determined from the MC samples. The efficiency functions e and background distributions V 
for X s £ + £~ event, dilepton background, K*£ ± h T background and ip background are obtained 
from large MC samples. Different efficiency functions are used for DS-I and DS-II since the 
acceptances are slightly different. A common background PDF used for both data sets with 
the averaged parameters. The Wilson coefficients in the PDF for X s £ + £~ are fixed at their 
SM values since this background's contribution is negligibly small. The K*h + h~ background 
shape VK*h+h~ is obtained from K*h + h~ data, applying momentum- and angular-dependent 
lepton ID fake rate corrections. 

The input parameters in the fit are q 2 , cos 9 and M bc . Only A 9 /A 7 and A 10 /A 7 are allowed 
to float in the fit. We use smoothed histograms for efficiency and dilepton background 
distributions, and unsmoothed histograms for other event categories. The background PDF 
shapes are found to agree with the the distributions in the data sideband well. 

Prior to the fit, we measure the background subtracted integrated asymmetry. To check 
the possible bias due to the event selection, we perform the null asymmetry test with the 
K£ + £~ control sample. We obtain 

A bk g -sub (jB ^ K+i+r) = 0.09 ± 0.14(stat.). (4) 

(5) 

Note that ApB g_sub is not corrected for the g 2 -dependent efficiency. The result for B — > 
K + £ + £~ is consistent with zero as expected (FigjSfa)). We thus confirm that the selection 
criteria do not produce any bias. On the other hand, a large asymmetry is observed for 
B -> K*£+£~ (FigH»,(c)): 

A^' mh (B -> K*£ + £~) = 0.56 ± 0.13(stat.). (6) 
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We fit the q 2 and cos 8 distribution for K*£ + £ with the PDF described in the previous 
section. Using A 7 = —0.330, the best fit Wilson coefficient ratios are 



A 9 /A 7 = -15.3 It*, 

A w /A 7 = 10.3-^: (7) 
while, with A 7 = +0.330, the best fit ratios are 

A 9 /A 7 = -16.3 If;?, 

A l0 /A 7 = 11.1 lf°. (8) 

Figure El shows the fit results projected on the background-subtracted forward-backward 
asymmetry distribution. 

We study the systematic error due to uncertainty in the physics parameters, finite q 2 
resolution, efficiencies for low momentum particles and the signal fraction. We change the 
A 7 value within the range allowed by the branching fraction of B — > X s ^ |26| to obtain 
the corresponding systematic error. The bottom quark mass m& is varied by ±0.2 GeV. 
The systematic uncertainty associated with the choise of form factor model is taken from 
the difference in fit results using the models of Ali et al. and Melikhov et al. j2?J. We 
ignore smearing in q 2 in the fits, although the actual resolution is about 0.7%; the effect 
is measured using a toy MC study. The efficiencies to find and identify low momentum 
particles are slightly better in MC events, so we simultaneously reduce the efficiency for 
pion with p < 0.3 GeV, electron having p < 0.7 GeV and muon with p < 1.0 GeV to obtain 
revised efficinecy functions for signal and background PDFs; we then refit the data with the 
new PDFs. The change in the fit results correspond to the systematic error associated with 
the identification efficinecy. The signal fraction is determined from a fit to M^ c distribution. 
We change the PDF shape parameters and take the difference as a systematic error: the 
parameters are modified by ±lo" for signal and dilepton background, by 100% for cross- feed 
events and ip background, and 10% (40%) for K*h + h~ background for muon (electron) modes. 
For the uncertainty in K^l^h^ background, we change the fraction of K*£ ± h T to dilepton 
backround, fxH^h^, by 20%. Table HP summarizes the contributions to the systematic error. 
The dominant contribution is from model dependence. 



source 


negative A 7 solution positive 


A 7 solution 




Ag/A 7 


A w /A 7 


Ag/A 7 


A w /A 7 


A 7 
m b 


+0.29 
-0.03 
+0.69 
-0.68 


+0.01 
-0.03 
+0.45 
-0.46 


+0.13 
-0.27 

±0.63 


+0.36 
-0.15 

±0.42 


model dependence ±0.66 


±1.72 


±1.04 


±2.23 


q 2 resolution 


±0.28 


±0.39 


±0.28 


±0.39 


efficiency 


±0.08 


±0.03 


±0.10 


±0.06 


signal fraction 


+0.43 
-0.47 


+0.22 
-0.33 


+0.43 
-0.46 


+0.37 
-0.40 


total 


+1.12 
-1.10 


+1.83 
-1.84 


+1.33 
-1.36 


+2.36 
-2.34 



TABLE II: Systematic error for each source. 
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FIG. 3: Forward- backward asymmetry distributions for K£ + £ and K*£ + £ . The points show the 
results obtained from M^ c fits while the curves show the fit results. 
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The best-fit results for the Wilson coefficient ratios for negative A 7 , 



and positive A 7 



A 9 /A 7 
A w /A 7 



Ag/A 7 

A 10 /A 7 



-15.3 tli± 1.1, 
10.3 1 5 3 2 5 ± 1-8; 



-16.3 ±f;? ±1-4, 
11.1 ±f;g ± 2.4, 



(9) 



(10) 



are consistent with the SM values Ag/A 7 = —12.3 and A\q/A 7 = 12.8. In Fig. 0J we draw 
contour lines in the Agj A 7 -A\gj A 7 plane based on the fit likelihood smeared by systematic 
error (assuming Gaussian). We also calculate an interval on AgAig/A 7 at 95% CL for any 
A 7 value, 



- 1401 < A 9 A W /A 2 7 < -26.4. 



(11) 



We determine the sign of AgAio to be negative, and exclude solutions in the first or third 
quadrant with more than 95% CL. Both second and fourth quadrant solutions are allowed, 
so we cannot determine the sign of A 7 A W yet. Figure El shows the comparison between the 
fit results for the negative A 7 solution projected on the forward-backward asymmetry and 
the forward-backward asymmetry distributions for several input parameters. We exclude 
new physics scenarios shown by the red and magenta curves, which have positive AgA w . 
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FIG. 4: Confidence level contours for negative and positive A 7 . Black, red, green, blue and black 
curves show 1 to 5 a contours. Black and red points show the best fit and the SM value. 

In summary, we have measured for the first time the ratios of Wilson coefficients in 
B — > K*£ + £~ by fitting to a normalized double differential decay width. We observe a 
large forward-backward asymmetry in B — > K*£ + £~. The fit result is consistent with the 
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FIG. 5: Fit results for the negative Aj solution projected to forward-backward asymmetry(solid 
blue) and forward-backward asymmetry curves with several input parameter including efficiency 
effect; A7A1Q sign flipped (dashed green), both ^7^10 and ^9^10 signs flipped (dash-dot red) and 
sign flipped(dotted magenta) to SM value. The new physics scenarios shown by the red and 
magenta curves are excluded. 

SM prediction and also with the case where the sign of A 7 A W nipped. We exclude positive 
A g A 10 solutions at more than 95% CL. 
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